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Abstract Explicit solutions of the two-dimensional floating body problem 
(bodies that can float in all positions) for relative density p ^ | and of the 
tire track problem (tire tracks of a bicycle, which do not allow to determine, 
which way the bicycle went) are given, which differ from circles. Starting point 
is the differential equation given by the author in [T2l [8] . The curves are also 
■ trajectories of charges in a perpendicular parabolic magnetic field. 

O 1 Introduction 

t> . 

In this paper we consider a class of curves, which are solution to three problems: 
*"^5 ' (i) the two-dimensional version of the floating body problem asked by Stanislaw 

O , Ulam in the Scottish BookpQ (problem 19): Is a sphere the only solid of uniform 

c/3 ' density which will float in water in any position? Thus we ask for the non- 

circular cross-section of a long cylindrical log, which can float in all orientations 
with the cross-section perpendicular to the surface of the fluid, 
(ii) It also deals with the tire track problem, which asks whether it is possible, 
to have tire tracks of the front and the rear wheel of a bicycle, which do not 
^> , allow to determine which way the bicycle went. 

j_i ■ (hi) Finally the boundary of the cross-section describes the trajectory of a charge 

moving in a perpendicular parabolic magnetic field. 

This paper is an improved and extended version of the solution given in [2] . 
It is extended insofar as figures for several curves are given. The discussion of 
the tire track problem is extended and explicit results for the cross-sections are 
given. Moreover the solution is formulated in terms of the four extreme radii, 
which makes the solution and symmetries more transparent. 
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(1) 



Two observations [3HH El [3 [S] are on the basis of the floating body problem: 

(i) Let h be the distance between the centers of gravity of the parts of the body 
above and below the water-line. Then the potential energy of the body is a 
constant times this height. Thus h has to be constant. 

(ii) If one rotates the body infinitesimally, then the line between the two centers 
of gravity will only remain perpendicular to the water-line, if the length 2£ of 
this water-line obeys 

3 \Ai A 2/ 

where At and A 2 are the cross-sections above and below the water-line, which 
of course are constant. Thus the length 2£ of the chord across the cross-section 
along the water-line is constant. Moreover since the areas above and below 
the water-line stay constant, the envelope 7 of the chords are touched in the 
middle of the water-lines. Since for relative density (density of body divided 
by the density of the fluid) p = 1/2, the water-line is the same, when the body 
is rotated by 180°, Auerbach|3j could give the solution for this case. Then the 
boundary of the cross-section is obtained by starting from a Zindler curve 7 
for the envelope, that is a curve which closes after the direction of the tangent 
to the curve has increased by n (such curves have typically an odd number of 
cusps, for examples see red curves and water-lines in figs. [TU1 to [T2l and [T5]) . If 
one adds tangents of length £ in both directions to all points of the curve 7, 
then these end points yield the boundary T of such a cross-section, provided 
this boundary is sufficiently convex. (By sufficiently convex it is meant, that all 
water-lines cross the boundary twice, but not more often) . In the following the 
case p ^ \ will be mainly considered. 

Another important property is obtained as follows: Consider two close-by 
water-lines A^A 2 and B^B 2 through the cross-section (Fig. [2 The body is kept 
fixed and the direction of the gravitational force is rotated.) Put the x-axis 
parallel to A\A 2 . Then the vector A\B\ is given by (daii, — £\d<$>) and A 2 B 2 
by (dx 2 , £ 2 d(j))- Constant length d£ = implies dx 2 = dx±. Constant areas 
At and A 2 imply d/ x = \i\d(f) = df 2 = \£\d(\>. Thus £t = £ 2 = £■ This 
implies, that the infinitesimal arcs at the perimeter dut = \J [dxt] 2 + £\{d(j)) 2 
and du 2 = {dx 2 ) 2 + £ 2 {d(f)) 2 are equal. Thus the part of the perimeter below 
the water-line is constant. Of course the same is also true for the part above 
the water-line. 




Fig. 1 Boundary and 
water-lines 



One can conclude the other way round: If a curve has the property that if one 
moves from two fixed points A\ and A 2 by constant arcs u along the perimeter 
and the length of the chord remains fixed, then also the areas separated by the 
chord stay constant. One argues that since £ stays constant one has dxi = dx 2 . 
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Since du% = du 2 also l\ — I2 and thus d/i = d/ 2 - In the main part of the 
paper this property will be used. By measuring the distance along the arc from 
a fixed point of the boundary one introduces the arc parameter u. It will be 
shown that for certain differences 2 8u of the arc parameter the length of the 
chord is constant, that is the distance between the point at arc parameter u— Su 
and at arc parameter u + Su does not depend on u. I call this the property of 
constant chord length. 

In the case of p = 1/2 the ends of the water-lines differ by an arc parameter 
which is one half of the circumference. Similarly one may look for solutions 
where the difference 2 5u of the arc parameter is a rational fraction of the 
circumference. Since however the area bounded by these water-lines and the 
part of the circumference is constant (the total cross-section is constant and 
the parts below the water also), the shape of such an equilateral triangle or 
quadrangle is constant and does not yield interesting results. An equilateral 
pentagon, however, of constant area and constant side length can vary its shape. 
Montejano er al.0 3] have used this property to obtain solutions for these 
special cases. Unfortunately there was no sufficiently convex boundary among 
the solutions. 

Another interesting property is that since l\ — £2, one has dyi — —dy2, 
which implies that the angles S between the tangents and the chord are equal, 

<*i = s 2 . 

The floating body problem is related to the tire track problem ( [SJ [HI fTU] 
and papers cited therein) . It goes back to a criticism of the discussion between 
Sherlock Holmes and Watson in The Adventure of the Priory School\Tl\ on which 
way a bicycle went whose tires' traces are observed. Let the distance between 
the front and the rear wheel of the bicycle be I. The end points of the tangent 
lines of length t to the trace of the rear wheel in the direction the bicycle went 
yields the points of the traces of the front wheel. Thus if the tangent lines in 
both directions end at the trace of the front wheel, it is open which way the 
bicycle went. Thus curves 7 for the rear wheel and T for the front wheel are 
solutions for such an ambiguous direction of the bicycle. The tire track problem 
consists in finding such curves T and 7 different from circles and straight lines. 
It is equivalent to the two-dimensional floating body problem with the exception 
that the tire tracks need not close or may wind several times around some point. 

Starting from r = ro(l + e cos(n^)) + 0(e 2 ) for the boundary of the floating 
body and carrying the expansion in e through to seventh order the author con- 
jectured that one obtains solutions for n — 2 different densities p for the same 
boundary [T2l 18] . This means that one and the same boundary T has the above 
mentioned property for several different chord lengths 21. Assuming now that 
this is also true, if the curve does not close as it is allowed in the tire track 
problem, one may assume that it also holds if the curve closes apart from an 
infinitesimal small angle difference Sx with chords of infinitesimal length 21. 
Then a differential equation of third order for this curve can be derived. After 
one integration this differential equation yields 
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r 2 _|_ 2r' 2 — rr" 
(r2 + r , 2) 3/2 =4^ + 26, (2) 

in polar coordinates (r, ^£>) with r function of ip and the prime indicating differ- 
entiation with respect to tp. k is the curvature of the curve, a is given by the 
limit = lim^ 3 jS\ for vanishing I and 8\ (defined by % — x of sect. I6.ip . b 
is the integration constant for this first integral. Thus these curves describe the 
trajectories of particles of mass m, charge e and velocity w in a perpendicular 
magnetic field, which depends quadratically on r, 

, , raw o 
B{r) = — (4ar 2 + 26). (3) 

Another integral of the third-order differential equation is 

f -4- r" 2c 

-2ar+ — . (4) 



lj.2 _|_ r /2^3/2 

If we introduce the inverse radius r = 1 /f, then this equation can be rewritten 
_,_ , f 2 + 2f' 2 - ff" 2a „ 

K ( r )^ (f 2 +f ;2)3/2 =-^4 +2C - ( 5 ) 

Thus f = 1/r describes the tractory of a charge in a perpendicular magnetic 
field 

- , , raw ( 2a \ 

B(f) = — ^ - 2cJ . (6) 

Elimination of r" from eqs. (J2j [4]) yields our basic differential equation 

1 



,.2 



VT 
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(7) 



It will be shown that the resulting curves solve the floating body problem, 
provided they are sufficiently convex and closed. Indeed there are such solutions. 
For the tire track problem the restrictions are less rigorous. But among them 
there are solutions, which require artistic mastery of the bicyclist, since he/she 
has to move back and forth. In more difficult cases the bicycle has to allow the 
steerer to be rotated by more than 180° degrees. 

Quite generally the following remarkable property of these curves will be 
shown: Consider two copies of these curves (not necessarily closed). Choose an 
arbitrary point on each curve. Then there exists always an angle by which the 
two curves can be rotated against each other, so that the distance between the 
points on the two curves stays constant, if one moves from the given points on 
both curves by the same arc distance. I call this the property of constant dis- 
tance. This remarkable property made it possible to obtain the above-mentioned 
differential equation by choosing the two points infinitesimally close. 

In [12l [8] the following necessary condition was shown: If for some chord 
the angles 6 between the chord and the tangents on the curve obey Si = S2 at 
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r\ 7^ T2 then also the derivative d< " Sl d ~ S2 ^ vanishes. This is a necessary condition, 
but it is not sufficient. I am indebted to Serge Tabachnikov for pointing this out 
to me. Here the differential equation |(7J) will be solved explicitly as a function 
of the arc parameter u, and it will be shown that the curve has the desired 
property. 

The outline of the paper is as follows: In sect. [2] an appropriate parametriza- 
tion is introduced. Depending on the parameters a, 6, and c one may have no 
or one or two (real) branches. Two branches means that the curves are not 
congruent. The equations are given in terms of the arc parameter. In sect. [3] 
the radius r as a function of the arc parameter u is determined. It is given in 
terms of Weierstrass p-functions. The invariants gi and g% are real. In the case 
of one branch the discrimant is negative, otherwise positive. This is discussed 
in sect. 0] In the following section the polar angle and thus the whole curve 
is determined as function of u. Starting from the differential equation ([7]) the 
solution can be written 

z(u) :=x + iy = Ce* ^ ~ 3 ^ e a « a «* (8) 
a(u + v) 

where u is the (real) arc parameter. The purely imaginary v and the real 
invariants g 2 and g 3 of the Weierstrass integrals C and c, and the constant C 
depend on the coefficients a, b, and c. \ is an arbitrary angle which defines the 
orientation of the curve around the origin and can be chosen arbitrarily. 

Instead of solving ([7]) one could immediately start with the curve defined by 
ijHJ and first show that the arc-parameter u indeed measures the distance along 
the curve, which is done in subsect. 15.31 

Secondly the above mentioned property of constant distance is shown in sect. 
O Given a difference 5u , there exists a difference of angles, by which two copies 
of the curve described by z and z have to be rotated against each other, so that 
the distance It = |z(u+ 5u) — z(u— Su)\ does not depend on u (eq. I141|) . If the 
difference of the angles is such that the two curves fall onto each other, then one 
obtains the desired property of constant chord length for floating bodies and tire 
tracks. The property of constant distance holds even between the two different 
branches in the case of positive discriminant. A remark on the expression of 
Finn for the curvatures of the curves is added in subsection !6.3l Finally the area 
above and below the water-line is determined and it is shown that the centers 
of gravity of the parts of the cross-section above and below the water-line lie on 
circles (compare Auerbach[3]). 

In the following section the periodicity of the curves is discussed and figures 
of several curves are shown. In addition a number of these examples are shown 
as animations on the internet [1 3 . These examples are indicated by an asterisk 
at the figure caption. 

The limit case, where the discriminant vanishes, is considered in section[5] In 
this case the curves can be expressed in terms of trigonometric and exponential 
functions. Moreover there is a simple construction principle of these curves 
explained in subsection 18.41 In section [9] the limit is considered, in which the 
curve oscillates around a straight line. (This does not constitute a solution of 
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the floating body problem, but of the tire-track problem.) There one finds an 
infinite set of Su s, for which the length of the chord between the points of arc 
parameter u — Su and u + Su is independent of u and where 2 Su is not a period 
of the curve. The curves found in this section are also trajectories of electrons 
in a linearly varying magnetic field as considered by Evers et al|14|. 

In the following section the carousels by Bracho, Montejano and Oliveros 
[HI [4] are discussed, where five chords or line segments close to an equilateral 
pentagon. The solution given here allows to reconstruct these curves. Some of 
them will be shown and discussed. 

I cannot claim that the solutions given here are the only solutions to the 
floating body problem, but I am not aware of any other solutions. There remains 
the question whether starting from the known solutions small deformations can 
yield new solutions similarly as a class of solutions were found from deforming 
the circle [8] . The solutions of the tire track problem are not restricted to 
closed and not terminating curves. In view of the paper by Finn [9j [10] who 
describes a procedure to continue from an initial back-tire track segment one 
expects a large manifold of solutions for this latter problem. 



2 Basic Equations 
2.1 A first discussion 

The curves described by the differential equation ([7]) will now be considered quite 
generally. For this purpose the angle <f> is introduced indicating the direction of 
the tangent to the curve with respect to the radial vector by 



dr 



r tan < 



(9) 



as shown in fig. [2j Then eq. ([7]) can be written 




Fig. 2 Figure indicating the angle <f>, 
part of the curve (heavy line) and the 
tangent to the curve. 



aq 2 + bq + c = y^q cos (j), (10) 
where the square q = r 2 of the radius is introduced. Thus the curves obey 



y/q < aq 2 +bq + c< *Jq. 



(11) 
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Whenever 



aq 2 +bq + c = ±^/q 



(12) 



is reached, the curve has found an extreme distance from the origin. 

Depending on the constants a, 6, and c the inequality (jlip allows for zero, 
one or two curves. This can be seen by plotting the parabola C depicting ±^/q 
and the parabola \J showing aq 2 +bq + c versus q. To the extend segments of the 
parabola \J lie inside the parabola C one obtains curves. These segments are 
indicated by thick lines. Various possible cases are shown in fig. [3l Figure [3^ 




Fig. Hi Fig. \B> Fig. He 




Fig. \3jg Fig. \%i Fig. \3jj, 

Fig. 3 Various cases of the relative position of the two parabola 



does not allow for any real curves, since the parabola \J lies completely outside 
C- In fig. [3b the parabola \J touches <. Only one radius is allowed, which 
corresponds to a circle. Fig. [3fc allows for some variation of the radius. It yields 
a curve oscillating between the two extreme radii. From this case one obtains 
good solutions for the floating body and the tire track problem. In fig. [3ji 
the parabola V touches the g-axis. Thus the curve shows a tangent in radial 
direction. In fig. [3^ the curve goes backward in ^-direction close to the smallest 
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radius. Fig. |3J allows for a small circle plus a curve. Fig. (3^ is the limit 
case where the curve approaches asymptotically the circle. Fig. [3)i allows for 
two curves both of which approach asymptotically a common circle, one from 
outside, the other from inside. Fig. [3J is the generic case of two curves. 

The extreme distances from the origin are reached, when cos cf> — ±1. This 
yields 

arf + brf +c-n=0 (13) 

where for real ri the \ri \ are the extreme radii of the curves, r,; is chosen positive, 
if cos = +1, and negative, if cos = — 1. The polynomial (|13[) can be expressed 
in terms of these extreme radii r,, 

4 



ar 4 + br 2 + c - r = aT\(r - n). (14) 

i=l 

Comparing the coefficients of the polynomial one obtains 

= ]>>, (15) 

% 

b = «E^ = i(E^) 2 -lE^ 2 = -lE^ ( 16 ) 

-1 = -a(rir 2 r 3 + rir 2 r 4 + rir 3 r 4 + r 2 r 3 r 4 ) = -aP r , (17) 
c = ar\r2r 3 ri — aP, (18) 

which allows the coefficients a, 6, c to be expressed by the r^. Thus the sum of 
the four r» vanishes. Here 

- r? (19) 

is introduced and use is made of the polynomials P r and P defined in appendix 
[B]in eqs. (|357[) and ([365j) and of eq. (|359p . In particular one obtains 

o=^-. (20) 

In the following a is assumed to be positive. Reversing simultaneously the signs 
of a, &, and c corresponds to reversing all rj, which yields the same curves, but a 
different orientation. If a = 0, then according to |2]) the curve is a circle (6 ^ 0) 
or a straight line (b — 0). 



2.2 Parametrization according to the shape 

Multiplication of all by the same factor yields similar curves, but no new 
shape. Since the four obey the restriction r,; = 0, only two parameters are 
left, which determine the shape of the curves. Thus the parametrization 

r 4 , 3 = r (l ± e), n, 2 = r (-l ± e) (21) 

is introduced, where ro sets the scale of the curve, whereas e and e determine 
its shape, e and e can always be chosen either real or purely imaginary. A sign 
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change of e and/or e corresponds to a permutation of the r-i. Therefore the e 2 , e 2 
plane is plotted in figured] Both e 2 and e 2 are real, ro is assumed to be real. 

Consider the four quadrants in fig. |4] e is real and e = —ifj, is imaginary 
in the left upper quadrant A. Here one has one real curve. In the right upper 
quadrant B both e and e are real. There one obtains two branches. In the 
right lower quadrant C e is imaginary and e is real. By exchanging e and e one 
comes back to quadrant A. Finally in the lower left quadrant D both e and e 
are imaginary. This does not correspond to any real curve. Thus it is sufficient 
to consider the quadrants A and B. 

Without restriction of the general case 

e-2>e>0. (22) 

can be required in region B. This corresponds to the ordering 94 > 93 > <?2 > Qi 
and is covered by the region Bi. Fig. [5] repeats fig. [3J with the assignments of 
the Ti. 




Fig. 5 The r t in Fig. 0i. 



The parobola in figure [J] 

[(e - 2) 2 - e 2 ] [(e + 2) 2 - e 2 ] = (e 2 - e 2 ) 2 - 8(e 2 + e 2 - 2) = (23) 



9 



and the line e 2 = e 2 decomposes the whole quadrant Be 2 > 0, e 2 > into 
six regions. If one allows an arbitrary assignment of the to ro(l ± e) and 
J~o(— 1 ±e) then one obtains 24 equivalent permutations. Since only the squares 
of e and e are plotted, this reduces to six sets of (e 2 ,e 2 ), which are represented 
by the six regions Bi to Bq. 



2.3 Equations to be solved 

From equation ([7]) one obtains 

dr \ ' 



dip J (ar 4 + br 2 + c) 2 



r 



2 



(24) 



and thus 

<ty = aq 2 + bq + c 

d<7 2qy/q- (aq 2 + bq + c) 2 

with q = r 2 . If one denotes the distance on the curve from some fixed point 
(arc-parameter) by u then one obtains 

pL = ^2~Tr^ = = i - (26) 

dip ar z + o + cr z aq + b + cq 1 



and 

du 



dg 2^g- (a<? 2 + bq 
First (sect. 13]) q(w) will be calculated from 



(27) 



(^) 2 = - 4 « 2 ri(^)> (28) 
dq 



(29) 



since 

(aq 2 + bq + c) 2 — q = (ar 4 + br 2 + c — r)(ar 4 + br 2 + c + r) 

~*J[{r-n)J\(r + r i ) = a 2 J\( q -q i ), (30) 



a 



Secondly ip(u) (sect. [5]) is determined from 

ip(u) = J dufaq+b + cq' 1 ). (31) 

Using eqs. ([5J 116[) the curvature is given by 

K = Aar 2 + 2b = a(Ar 2 - ^ r 2 ). (32) 
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One concludes that a convex curve oscillating between ro(l + e) and ro(l — e) 
has to fulfill 

4e<e 2 -e 2 . (33) 
This condition guarantees for positive e that the curvature does not change sign. 

3 Radius as Function of Arc Parameter 

The solution of the differential equation (|28[) reads 

9 = 94—-^ (34) 

P{u) - Pi 

with the Weierstrass p-function. To show this, one calculates the derivative 



dq p'(u) 
— = -qi(Pi -Paji—r-, v 

du (p(u)-pi) 

2^4 (Pi -P3) 



94(Pi-P3) 7r7 ^TT2 (35) 



, . Y2 V (P( u ) ~ e i)(p( M ) ~ e 2)(p(u) - e 3 ). (36) 

From eq. (fM)) one obtains 

1 _ 1 /■? _ A p(w) - Ci _ pi - e,; ( q p 3 - e,- 



p(u)-pi Pi-P3\94 / p(u)-pi Pi-P3\94 Pi-ej 

(37) 

Eq. is fulfilled provided the relations 



A 2 4 n i= l(Pl- e i) P3-Ci ■ , A , Q ca 

4a z = L \\ r^— , ft =94 , «7^4 (38) 



hold. The last eqs. and Y^=i e i = yield four linear homogeneous eqs. for pi, 
Ps and the ej. Together with the first equation one obtains the solutions 

a 2 

pi = a 2 <? 4 (-<74 + 9i + 92 + 93) - irPra, (39) 
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a 2 



2/ 919233 . . . \ " D IAC\\ 

P3 = a{ 1- 9i92 + 9i93 + 9293) - -5-fVn, (40) 

94 o 

a 2 

ei = a 2 (9i94 + 9293) - y-Pm, (41) 

a 2 

e 2 = a 2 (g 2 94 + 9i9s) - yP m , (42) 

a 2 

C3 = a 2 (9394 + 9i92) - y-Pm (43) 

with the polynomial P m as defined in eq. (|358|) . Permutations of the r*j yield 
permutations of the ej. Therefore the invariants <?2 = — 4(eie2 + eie^ + e^e?,) 
and 53 = 4eie2e3 and also p(w) are invariant under all permutations of the r,. 
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For future purposes the parameter v is introduced by 

Pi = p(v), (44) 

which yields 

p'(v) = T2ia 3 {q 4 - qi )(q 4 - q 2 )(q 4 - q 3 ), (45) 
p (v) = 2a A (q 4 - q 4 ){q 4 - q 2 ){q 4 - qi)(?>q 4 - qi- q 2 - q?,)- (46) 

From these expressions one concludes 

a 2 a 2 

p{2v) = -—{ql + ql + qi + qD + ^P™, (47) 

p'(2v) = T 2ia 3 P q =T^- = T2ia (48) 

with the polynomial P q defined in (|356p and by use of eqs. ([2171 I36ip . These 
expressions are invariant under all permutations of the r^. 
Using the condition ^\ = one obtains 

P3 = p(3v). (49) 

To show this one starts from 



o'{v) - p\2v)) 2 

i( P (v) - p(2v)y 



p(3v) = -P{v)-P{2v) + Z^, in a . (50) 



P{ v) = - P ( W )- P ( 2 ,) + yy + ^f)) a , (5i) 

4(p(u) - p(2u)) 2 



from which one concludes 



(p(») - p(2t>)) 2 



K&0 = P(«) - ^TTCT - (52) 



Comparison shows 

p(t>) - p(2t>) = a 2 ( - 2P + q 4 { qi +q 2 + q 3 - q 4 )). (53) 

with the polynomial P defined in eq. (|365p . By means of eqs. (|362l I366[) one 
finds 

p(v) - p(2v) = -2a 2 r 4 P r (54) 

and 

(p(v) - p(2v)) 2 = Aa\ 4 P q . (55) 

Insertion into eq. ([52)1 yields eq. ([49)) . One obtains the derivative p'(3v) by 
means of (|346|) 

p'(Sv) = -^p-p'(v), (56) 
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which by use of eq. (|366p can be written 



p'(3v) 

The solution (f3~4"| can be written 

/ 3u,-3w 
<?(«) = q±£ u 



P 

4 



p'(v). 



94" 



ct(m + 3v)a(u — 3w) t 2 (w) 
cr(u + v)cr(u — v) a 2 (3v)' 



(57) 



(58) 



where the arc parameter u is measured from the extreme radius r = r±. The 

notation £ n j" 1 j J for double-periodic functions of it is explained in appendix 

IA.21 Apart from the origin it has zeroes at Uj and poles at Vi and behaves like 
u n at the origin. Correspondingly the radius reads 



r(u) = C 



a(u + 3v)a(3v — u)\ 



a{u + v)a(v — u) 
The prefactor C is rewritten by means of the identity 

a(3v) 



C = r 4 



cr(3«) 



a(v)cr 2 (2i;) 



p(v) - p(2u), 



(59) 



(60) 



which can be seen by realizing that both the left and right hand-side of this 
equation is double-periodic and can be written 



[v 

One obtains 
C = 



2m i — 2oji 2aJ2 — 2aJ2 2a;i+2aJ2 — 2a;i — 2aJ2 2oji — 2oJ2 2u2 — 2oji 
3 ' 3 ' 3 ' 3 ' 3 ' 3 ' 3 ' 3 

U)\, — OJl, U)2, — W2j Wl + CJ2, — Wl — k>2 



(61) 



(62) 



where eqs. (1481 13441) have been used. Thus the factor r 4 drops out and the 
prefactor C is invariant against any permutations of the r^. 



~Pr 



\2v) 



a 2 {2v){p{v) - p{2v)) 2a 2 {2v) p'{2v)a 2 (2v) a(Av) 



4 Single Branches and Pairs of Branches 

In the preceding section the radius square q has been obtained as a function of 
the arc parameter u in eq. (|34|) . At u = the function p(u) diverges and q 
approaches q 4 — r 2 . As u is increased up to w 3 one obtains p(w 3 ) = e 3 and q 
approaches its minimum q = 94(63 — p§)/{&3 — Pi) — 93 along the curve. Then 
q increases until at u = 2^3 it returns to q = q 4 . Instead starting at u = with 
q = q 4 one can start at u = w 3 with q = q 3 . Therefore consider 

q{u + Wi) = q 4 — — ■ r -— , (63) 

p{u + uj t ) ~ p(v) 
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with half-periods UJi, i — 1,2,3. uj is called a half-period, if it is not a period of 
p, but 2uj is a period. The expression (|63|) can be rewritten 



q(u + uJi) — q^So [u + uj. 



94- 



3u, — 3u 
u, — w 

3w — u>i, ~3v — u>i 
v - UJi, -v - UJi 
e-i - P3 p{u) - p(3v - UJi) 



p(uji)-p(3v) . 

= ^~T~\ TT Ea u 

PK^t) - p{v) 



- Pi P{u) ~ P(V - UJi) 



p{u) -p\' 



with 



pf=p{v^), = p(3w«), (65) 
Explicit calculation by means of 

e? + e,-e fc + ejp(v) 

p(u-cJi) = ^— , (66) 

PW - e t 

where j and k are the elements out of (1,2,3) different from i, shows that Pi 

(i) 

and p 3 are obtained from p\ and p3 by exchanging and qn in the expressions 



([391) and Pjl. 

The expression f|47[) for p(2t>) is invariant under permutations of the 
Solving for v from this equation yields besides v from ((551) a l so the yW, i = 
1,2,3. In the following v will be chosen, so that Syp'(u) < 0. Then the upper 
sign in eqs. ()45l |48j) apply. Since from ()346|) one obtains with u = — cjj 

P ^ w )=- p( 7r ) ~ e< p , (^ w 

p{v) - et 

one finds that p'(y^) is obtained from p'(v) by exchanging g.; and 54 and one 
has also to take the upper sign. 

4.1 Single branch 

If two of the r.i are complex and two are real, then one (real) curve is obtained. 
The choice of complex r 2 — r\ and of real and r% yields 

92 = q*i, 94 > 93- (68) 
Then e 3 is real, whereas e\ and e 2 constitute a conjugate complex pair 

e-2 = ej, e 3 = 63. (69) 
Consequently the discriminant 

V = 16(e 3 - ei) 2 (e 3 - e 2 ) 2 (e 2 - ei) 2 < (70) 
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Fig. 6 Two elementary cells of periodicity for negative discriminant inside the 
rectangle. 

is negative. For negative discriminant T> the periods and half-periods are located 
in the complex plane as shown in fig. [S] 

The rectangle < < 2w 3 , < 3z < 2o//i covers two elementary cells 
for the double-periodic function p(z). The function is real along the straight 
lines drawn. For constant (horizontal lines) it varies between +oo and e%. 
For constant 5ftz (vertical lines) it varies between — oo and e 3 . The singularities 
are at the points indicated by •. p(») = oo. The derivative of the function 
p(z) vanishes at half-periods ui. There are three different half-periods which 
do not differ by periods. They are indicated by O , -J- , and □ . Thus one 
has jp'(o) = p'( + ) = p'(n) = 0. At these points the function approaches 
p(u>i) — e». Examples are p(u>i) = p{— lo\) = p(o) = e\ with 2u)\ = u>' — u> 3 , 
p(uj 2 ) = p(-oJi) = £>(+) = e 2 with 2uj 2 = -u' - u 3 , p(uj 3 ) = p(u)') = p(u) = 
e 3 . 

Next consider the location of p\. In section [3] the pi related to q± was 
considered. To indicate this connection one may write p\ = p^ and generally 
introduce 




(71) 



Here one is interested in the two real p\ and p\ . One obtains 

P^ - e 3 = -a 2 (q 4 - <Zi)(q4 - <?2) < 0, 
= -a 2 (<73 - <7i)(<?3 - 92) < 0. 



(72) 
(73) 
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Thus one has 

P^\pf<e s . (74) 

The real part of given by 

p(««)=pW (75) 

lies on the vertical lines in fig. [6] 

3?w (4) ,5Rw (3) =7iw 3 (76) 

with integer n. In the following the main choice for z/ 4 ) = v will be the one 
with 

5Ru = 0, < 3w < w', (77) 

so that 

3?p» = 0, 9p» < 0. (78) 

in agreement with the requirement introduced before. As mentioned before 
adding u 3 to u is equivalent to exchanging q 4 and q 3 . Both choices yield the 
same curve. To exchange q<± with q\ or qi does not yield a positive q and thus 
a real curve. 

The are obtained from 

^3 = 7= CJ =1 = , (79) 



'Ha a/ Ha 



#3 = V(e3-ei)(e 3 -e 2 ), (80) 

where for the elliptic integral K(m) the notation of Abramowitz and Stegun[T5] 
is used. The reader is refered to chapter 17 of this handbook for elliptic integrals 
and chapter 18 for Weierstrass elliptic functions (watch permutations of indices 
i of 6j and u>i in comparison to the use here), and to appendix |A1 of the present 
paper. 

4.2 Pair of branches 

If all four Ti are real, then one obtains two branches. The should be ordered 
according to 

<74 > <73 > 92 > q%- (81) 
Then the differences of the ei, e2, e 3 obey 

e 3 - e 2 = a 2 (q 3 - q 2 )(<74 - Qi) > 0, (82) 

e 2 - ei = a 2 (g 4 - <73)(?2 - 9i) > (83) 

and thus 

e 3 > e 2 > ei. (84) 
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Fig. 7 One elementary cell of periodicity inside the upper right rectangle for 
positive discriminant. 



In this case, where all three are real, the discriminant V, eq. (|70p . is 
positive. Then the periods and half-periods are located in the complex plane as 
shown in fig. [7] 



One elementary cell of periodicity is given by the rectangle < 5Rz < 2uj 3 , 
< 3z < 2u>i/i. Again the function is real along the straight lines drawn. The 
function approaches infinity at p(») = oo. The half-periods are indicated by 
□ , O and -f. The derivatives vanish at these half-periods, p'(n) = p'(o) = 
p'(+) = 0. The values of p at the half-periods are denoted by p(ua) = ch, 
p(wi) = p(o) = ei, p(w 2 ) = pO ) = e 2 , p(w 3 ) = p( □ ) = e 3 . 

One obtains the inequalities for the defined in (|7l|) 



M 4) 

(3) 
Pi 

(3) 
Pi 

(2) 
Pi 

(2) 
Pi 



ei = -a 2 (q 4 - 92)(94 - 93) < 0, 

ei = -a 2 (<72 - gi)(g3 - 9i) < 0, 

P2 = a 2 (q 3 - <? 2 )(<7 4 - 93) > 0, 

e3 = -a 2 (g 3 - <7i)(<?3 - 92) < 0, 

e 2 = a 2 (g 2 - qi)(q 3 - 9i) > 0, 

e3 = -a 2 (<73 - <?2)(<74 - 92) < 0. 



Thus one has 



(4) (1) „ ^ . (3) (2) 

pi < ex < e 2 < p\ \p\ 



< e 3 . 



(85) 
(86) 
(87) 
(88) 
(89) 
(90) 

(91) 
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Consequently one finds for the real part of defined in (f75j) 



5Rw (4) , %tv {1) = 2nuj 3 , Kw (3) , $tv {2) = (2n + l)w 3 (92) 

with integer n. In the following the main choice for = v will be the one 
with 

9fJv = 0, 0<3w<wi, (93) 

so that as in the case of one curve eq. ([75]) and the upper signs in eqs. ([^51 
I48p apply. Adding L03 to it again yields the same curve with u running through 
• and □ and the radius oscillating between ^fql and y/qH. However, adding 
u>i or UJ2 yields a different variation of the radius between ^Jql and ^fq?- It 
corresponds to u running along a horizontal line in fig. [7] through o and +, 
where p is again real. This is equivalent to subtracting u>i or L02 from v. Thus 
in this case one obtains two branches. 
The u>i are given by 

^3= , e3 ~ ei \ wi=i ; e3 ~ ei ' . (94) 
V e 3 - e i V e 3 - ei 

There are the special cases where two of the radii r.i coincide. In these 
cases the discriminant vanishes. In fig. [3f one has t\ = ra. This case yields a 
small circle plus a curve with varying distance from the origin. In this case the 
imaginary period goes to infinity. In fig. [5h, where r 2 = r% there are two curves 
asymptotically approaching the same circle. There the real period approaches 
infinity. The limit case is shown in fig. [3g, where r% — r 2 = r 3 . In this case 
both periods approach infinity. These limit cases are considered in section [5] 



5 Angle and Curve as Function of the Arc Pa- 
rameter 

5.1 The angle 

The angle i\) as function of u is determined from eq. (|31| 



*(«) = / du(a q + b + cq -i) = [du (a q4 ^^ + b + ^^t^l) . (95 ) 



With 



one obtains 



du . . 1 . . = ( 2«C(») + In f ) 1 I > 

p(u)-p(w) p'(v) \ \<r(u + v)J 1 



, a(v — u)\ , f a(3v — u)\ 
if}(u) =x + c u + c 1 ln[ + c 3 In (97) 



a(v + u)J \a(3v + u) 
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with the real integration constant x an d 



c = aq 4 + b+ — + 2ciC(V) + 2c 3 C(3v), (98) 
c a = ^rf, (99) 
C3 = ^7^- (100) 



CX=C3 = -i (101) 



Evaluation of ci and C3 yields 



where the sign convention (|78p is applied. Further one obtains from |T 

aq 4 + b+ — = — . (102) 

94 r 4 



^i(u) = x + c u - - In — — ■ — - - - In — — ■ — - (103) 



This yields 

2 \cr(u + u)y 2 Vo-(3?; + w) 

with 

co = --i(C(«) + C(3«)). (104) 
r 4 



Using 



one obtains 



C(2v±v) = a2v)±av) + 2 ^^ ) (105) 

2(p(2w) - p(v)) 



«*,) + «.) = 2C(2») + ^^l. _ 2C(2 „) _ . 2C(2 „) _ 1 (106) 

and thus 

c = -2iC(2«), (107) 

which yields 

\a(v + u)a(3v + u) J 

5.2 The curve 

Combining the result for the radius and the angle one obtains the curve as a 
function of the arc parameter u 

«.) := * + * - r(«)e*> = (109) 
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as claimed in flSJ). 

Changing simultaneously the sign of u and v does not change z, 

z(u,—v) = z(—u,v). (HO) 

Thus changing the sign of v is equivalent to changing the orientation of mea- 
suring the arc-parameter. 

Taking the complex conjugate of z yields 

[e~ ix z(u, v))* = e~ ix z{u, -v) = e~ ix z{-u, v). (Ill) 

Thus the curve has a mirror axis which runs through the origin and e lx . 

By means of eqs. (|3291 13301) one finds that adding a period 2w to v does not 
alter z, 

z(u,v + 2w) = z(u,v). (112) 

One obtains 

e i«/)(M+2ai 3 ) _ gi-i/ , ('") e i'0p<!r (113) 
gi^por _ e -8-uC("3)+4w 3 C(2«) (114) 

Thus the curve is periodic under rotations by -0p CT - This equation determines 
the angle tppcr only modulo 2n. Define 

V> c := 4tt + 8i;iCM - 4iw 3 C(2i;). (115) 

Then e^p" = e"^ c . In sect. [7] it will be shown that this angle equals ifj pel in the 
region A' of fig. There it will be discussed, how ^ per and ifjc are related in 
the other regions. 

Adding a half period to v yields 

e iV'pcr(f + w) _ e i-0pcr(>) e -8ia<(w3)+8iW3CM _ e i'0pcr(f) (H6) 

where use has been made of the Legendre- relation (|334|) . Thus two curves with 
v differing by a half-period have the same e 1 ^" 1 . 

5.3 More on the curve 

Instead of verifying that (| 1 09[) obeys the differential equation ([7]) one may im- 
mediately start with the expression 

z = x + iy = C ^±^)^u (n7) 

a(u + v) 

and demonstrate that it is solution to the problem by first showing that u is 
really the arc parameter provided one has the appropriate relations between v, 
v' ', and Co, which requires 

dz 



du 



1 (H8) 
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and secondly show that it is possible to find differences 2 Su of the arc parameter, 
so that the distance \z(u + Su) — z(u — 5u)\ does not depend on u. Here it 
will be shown that (|118[) is fulfilled with appropriate choices of v', Co, and C for 
given real g 2 , <73 and imaginary v. The investigation of the property of constant 
distance is postponed to sect. [5] One obtains 



d£ 

du 

dz* 

du 

for real u. Thus 



= z$(u), $(u) = C(u + v') - C(u + v) +ic , (119) 
z*$*(u), $*(u) = £{u + v'*)-C(u + v*)+ic* (120) 



dz dz* , , , . 

— = $ (u)zz* = 1 

du du 



has to hold for real u with 



= CC 



, o{u + v')o{u + v'*) cl{co _ c » )u 
a(u + v)a(u + v*) 



(121) 



(122) 



Since in this representation $(u)$*(u)zz* is holomorphic in u, it has to be 
constant for complex u. Since <!>(u) and are periodic functions in u, also 

zz* has to have this property, which yields 



$(u) = c u £ ( u 



Ui,U 2 

-v. —v' 



$*(u)=c* a £ [u 
= c z £ ( « 



1 /* 

— w , — V 



(123) 



where ui, U2 depend on the choice of cq, and c u and c z are constant prefactors. 
Further v' + v'* can differ from v + v* only by an integer multiple of 2^3. The 
product of the three functions in Q123P yields 



£ u 



ui, u 2 , u 1; u 2 
-v, —v, — v*, —v* 



(124) 



In order that this function is constant one chooses modulo periods 2oj 

Ul =u 2 = -v*. (125) 
(The choice u\ — ~v or u 2 = — v would yield infinite c u .) Thus 



$(u) oc £ I u 



-V ,-V 
-v, —v' 



which requires 
Thus one concludes 



-v + 2v* + 2u. 



-3%v + 23w. 



(126) 
(127) 
(128) 
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The real part of v can be absorbed into u. Thus one may choose v purely 
imaginary. Since the functions S are invariant against addition of a period to 
any of its arguments, one may choose u> = 0, which yields v' = — 3v. Since <£>(u) 
has to vanish for u = —v*, one obtains 

ic = 2C(2u). (129) 

can be written 

$(«) - C(« - 3«) - C(« + «) + 2C(2«) = - f^ 2 ^ r, (130) 

where eq. (|333j) has been used. The prefactors in eqs. (| 1 23|) are 

C = - C (3i,) + 2C(2i,)-C(«) = ^^, (131) 



Thus one has the condition 

a 4 (2v) 



c u c* u c z = -C 2 4S4 = 1 (132) 



in agreement with eq. J62 



6 Line Segments of Constant Length 

6.1 Property of constant distance 

In this section the distance 2£ between points at u + Su and u — 8u on the 
same or on different curves of equal or different branch is considered 

{2lf = (z*(u+ Su)-z*(u- Su)){z(u+ 5u)-z(u- Su)) (133) 

in order to find line segments of constant length 21 independent of u, 

Z ( U ) = CV* ^" 3 ^ 2 "^, (134) 
a(u + v) 

|( U ) = cj* a ^~ 3 % 2 "« 2 *\ (135) 
a(u + v) 

For simplicities sake assume that v and v are purely imaginary. Then the 
conjugate complex of z and z is obtained simply by reversing the sign of v, v, \ 
and X. Since u and Su are real, they are left unchanged upon determining the 
conjugate complex of z and z, and only after this conjugation they are extended 
into the complex plane. Thus one has 

Z*(U) = Ce -ix g (^ + 3 ^ e -2<(2.) ! (136) 

a(u — v) 

z*(u) = Ce-«^±M e -2«C(2*). (137) 
a(u — v) 
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The expression (21) 2 becomes an analytic function in u and Su . Choosing 

V = v + vu)\ (138) 

allows not only to consider line segments between copies of one and the same 
curve (even v), but also to connect two curves of different branches (odd v). 

Similar to eq. (|114[) addition of 2ui to the argument of u multiplies both 
z and z by the same factor e 8 «CM-4w£(2i>) an( j divides both z* and z* by this 
factor. Thus A£ 2 is a periodic function in u. 

Apparently z*(u + Su ) has a pole at u = v — Su and z*(u — Su) a pole at 
u = i+ Su . In order that (2£) 2 is constant, the other factor z(u+ 8u) — z(u— Su ) 
has to vanish for these values of u, 

z(v) — z(v — 2 Su ) = 

_ Ce ix e 2va2v) _ g e ix <~ 2 SU+V- 30) e _ 4 Su C(2»)+2, C (2*) = q (13g) 

a(-2Su + v + v) x ' 

z(v + 2 Su ) — z(v) — 

Cc iX g ( 2 6u + V ~ 3i; ) c 4 fa C(2t>)+2tiC(2t-) + ^ e ix e 2«C(2«) = Q ( 14Q ) 

cr(2 i5w + t> + i) 

Both eqs. yield 

e i(X"X) = ( !)l-^ c -2fa ((T(2t»HC(2C)) °~( 2 Su+V + V) ^ 

cr(2 foi — u — £>) 

Similarly z(u + Su) has a pole at u = — v — Su and z(u — Su) has a pole at 
u = — v + Su . It turns out that z*(u + Su) — z*(u — Su) also vanish, if (|141| 
is fulfilled. Thus if condition (|14ip between x~ X an d <5w is fulfilled, then the 
distance 21 does not depend on the arc parameter u. This proves the property 
of constant distance. 

Next the number of solutions will be determined. First one realizes that 
X — X is a monotonic increasing function of Su , which can be seen as follows. 
The derivative 



d(x ~ X) 



A Su 

can be rewritten 



2i(C(2u) + C(2v) - C(2 Su + v + v) + C(2 Su - v - «)) (142) 



d(x - X) = f ; P (2fa)-p(2») for 5 = U > 

H A?/ 1 -2ip / (2-»+a;i) p(2fa)-ei e f _ , 



(143) 



One convinces oneself that all numerators and denominators are positive. Thus 
X — X is a monotonic increasing function of <5u . Next consider the change of 
X — X during one period, that is while increasing Su by u>s, 

Ax = Aa,iga(25u +v + v) - Aargcr(2<5u -v-v) + 2iuj 3 (((2v)+((2v)). (144) 
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Since < 3(2u) < u'/2 and < 9f(2u) < w x one obtains from eqs. [[55111552]) 
A X - -2tt - 4i(« + fi)C(w 3 ) + 2i^ 3 (C(2^) + C(2t))) = 2tt - Mlt^M , (145) 



If curves are considered, which close after n periods, that is for which 

772 

ip c (v) = 2ir— 
n 

with m, n coprime holds, then x ~ X increases by 

n — m for v = v. 



nAx = 2ir x 



2n — m for v = V + U>i, 



since from eqs. (|115l) and (|334[) one obtains 

ip c (v + uoi) = ^c(«) - 47T. 



(146) 



(147) 



(148) 



The numbers n — m and 2n — m are the number of solutions Su for given x — X 
after they repeat, since with 5u also 8u + 71.CJ3 is a solution, which yields line 
segments between the same points of curves. If chords on one and the same 
curve are considered (v = v, x — x); which constitutes the original problem, 
then one solution is the trivial solution Su — 0. Thus one is left with only 
n — m — 1 non-trivial solutions, which moreover appear pairwise, since with Su 
also — Su is a solution. For m/n = 1/2, for example one has only the trivial 
solution. If however two copies of the curve are rotated against each other, 
then a set of line segments of constant length is obtained. Examples are shown 
in figures [5] and [5] The asterisk at the figure number indicates that there is 
a corresponding animation on the internet |13j. where one can see the chord or 
segment of constant length moving. 





Fig.* 8 m/n = 1/2, e = 0.2 



Fig.* 9 m/n = 1/2, e = 0.2 



Here and in the following the curves T are shown in black and blue, the 
chords and line-segments are shown in green and cyan. The color switches in 
the middle, where the lines touch the red envelopes 7. 
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6.2 Length of the line segments 

Next a simple expression for the length of the line segments is derived assuming 
that Su has been determined. Since 4£ 2 is independent of u, one considers the 
expression in the vicinity of diverging z*{u + Su ), that is around u = — Su + v. 
There the residuum is given by 

res = Ce-' lx a{Av)e- 2va2v) . (149) 

This divergence is compensated by the zero of z(u + Su ) — z{u — Su). In order 
to obtain 4£ 2 the residuum has to be multiplied by the derivate 



dz(u + Su ) 



du 



dz(u — Su ) 



«=- Su +v du 



(150) 



u— — 5u +u 

= $(v)z(v) - §>(-2Su +v)z(-2 5u +v) (151) 
Using (|130p one obtains 

*(«)=0, H-28u+v)= ~fft\ FV ( 152 ) 

p[2v) — p(—2 Su + v — v) 

Thus the residuum has to be multiplied by —$(—2 Su + v)z{—25u + v), and 
since z(— 2Su + v) = z(v), one obtains 

At 2 =res-(-z(v)${-2Su +v)) - -C 2 a{4v)—— -. (153) 

v w v " v ' p{2v) - p{-2 Su + v - v) v ; 

Making use of eq. ([B^)) and that 2v and 2v differ by a period 2u>, one obtains 
finally 

i 

U 2 = (1541 

p(2 <5u + « - v) - p(2v) ' v ' 

This expression yields bounds on the length of the line segment. Within one 
and the same curve one uses 

e 3 < p{2 5u) < oo (155) 

For p(2 Su ) — oo one has 4f 2 = 0. In the other limit case one finds 

a 2 11 
e 3 - p(2v) = —(q 4 + q 3 -q 2 - qi) 2 = -, ; = -, ; , (156) 

where eqs. ([2TJ)) , p57p . and p60p have been used. This yields the bounds 

Q<2t< |n + r 2 | =r 3 + r 4 . (157) 
In the case of two branches being connected by a line segment one has 

d < p(2Su +lox) <e 2 . (158) 
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Then the bounds read 

|ri + r 3 | = r 2 + r 4 < 21 < |r 2 + r 3 \ = r x + r 4 . (159) 

There are obvious geometric bounds. For line segments between points on 
the same curve one has the upper bounds r 4 + |r*3 1 and if there is a second branch 
also |t~2 I + In. |- Lower bound is 0. For line segments between different branches 
one has upper bounds r 4 + |rj,| and |r*a | + |ra| and lower bounds r 4 — |ra| and 
| r~3 j — ] r*i | . It turns out that the limits given by eq. (|154[) are in some cases more 
restrictive. They are given in the tables (|160|) and (|161|) . 
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6.3 Remark on the approach by Finn 

David Finn considers the bicycle problem in [HI [TO] ■ Denoting the curvature of 
trace of the rear wheel by up he obtains for the curvature of the trace of the 
front wheel 

^ ±Su) = (iHZ)^ ±e rS^> (162) 

where the two signs apply to the bicycle riding in both directions. He gives 
ftp for an intervall 2 Su and continues so that this equation is fulfilled. If one 
eliminates up then one obtains 

d(n(u + Su ) + k(u — Su )) . . . . , . 

-!— i L i LL = (k(u+5u)-k(u-Su)) 

du 



^ - J(k(m+ Su) +k(u- 5tt)) 2 (163) 



Performing some algebra one sees that this equation is fulfilled by our solution. 
To do this one may express n and i in terms of p(u), p( Su ) and p(v). 



26 



6.4 The area 

Next the area A 2 below the water-line is determined, which is the area of the 
sector ^4 S from the origin to the end points of the water-line with arc parameters 
u — Su and u + Su minus the area At of the triangle with corners at the origin 
and at the end-points of the water-line, 

A 2 =A S -A t . (164) 

The area of the sector is obtained by integrating 

d.4 s = — (xdy — ydx) = — (z*dz — zdz*). (165) 



Use of dz — z$(w)du, dz* = z*$*(u)du, zz*$(u)$*(m) = 1 from section IQl 
yields 

dA s = i***(*(u) - **(u))du = i - du. (166) 

Insertion of eq. (|130p and use of eq. (J35J) yields 

A s = / — (2p(2«)-p(« + «)-p(«-i;)) (167) 

Ju- Su Sa 

= ^-U8up{2v) + £{u + Su + v) + ((u + Su -v) 

8(2 

-£(u- Su +v)-((u- Su -v)) (168) 

= — (4<5up(2w) + 2C(5u +v) + 2((Su - v) 
8a 

p'(Su+v) p'(Su-v) 



(169) 

Su + v) — p(u) p( Su — v) — p(u) / 
The area A peY of the sector by increasing u by 2uj% , that is Su = is given by 

Aper = ^( W 3P(2«) + CM). ( 17 °) 

Then the density p for a cross-section with ^ per = 2-7r/n is given by 

A 2 



(171) 



The area of the triangle spanned by the origin and the points on the curve at 
arc parameters u + Su and u — Su is given by 

At = -(^(w — Su)y(u + Su) — x(u + 5u)y(u — Su)) (172) 

= — - {z*{u — Su )z(u + Su ) — z*(u + Su )z(u — Su )) (173) 

C 2 ( o{u — Su +3v)a(u+ Su — 3v) 4j u £pj„) 
4i \ a[u — Su — v)a(u + Su + v) 
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a(u + Su + 3v)a(u- Su - 3v) c _ 45uC(2v) \ 
a{u + Su — v)<r(u — Su + v) J 
CP_ f a 2 (5u - 3v) p(u) - p(5u - 3v) c4SuC(2v) 



4i \ a 2 { Su + v) p{u) — p(5u + v) 



a 2 (Su + 3v) p(u) - p(Su + 3v) c _ 4SuC(2v) \ 
a 2 (5u — v) p(u) — p(Su—v) 



In the last expression one rewrites 

a 2 (5uT3v) <j(28u^2v) p'( Su ± v)a(^f4:v) 



a 2 (5u ± v) cr(2 Su ± 2w) p{ Su =F 3u) — p( Su ±v) 



(176) 



where numerator and denominator on the left hand side have been multiplied 
by a 2 (Su ± v) and eqs. ((3271 l344f have been used. Use of eqs. (|i8 l l62 ]) yields 

A = A + P'( Su + v ) n + n+ p'(Su-v) 
2 8a(p(Su +v)-p(u)y J + ) 8a(p(Su -v)~p(u)y J ~> 

(177) 

with 

A = —{A8up{2v)+2Q(5u+v)+2a8u-v) (178) 
8a 

p'{5u +v) p'{Su -v) 

+J+~n n — ; — \ + J- 



p{ Su — 3v) — p{5u + v) p( Su + 3v) — p{8u — v) 
_ <j(2Su T 2v) ±4SuC(2v) 
J± ~ a(2Su±2v) ■ 1 ' 

In order that the area A 2 does not depend on u, the factors 1 + f± have to 
vanish. This yields exactly the condition (I141j) with X = Xi v = Q, v — v for 
the chord of constant length. This was expected. The constant area below the 
water-line is then given by A 2 = Aq with f± = —1. The expression for the area 
„4o can be simplified. One obtains 

" {<) " ±r] =2((8u±v)-({2 8uT2v)T(('iv) (180) 



p( Su =F 3v) — p(Su ± v) 

by using theorem I of appendix lA.21 The functions on both sides of the equation 
are periodic in Su , their singularities within an elementary cell are simple poles 
at Su — =pw with residua 2, and at Su — ±v and Su = ±v + Ui with residua 
— 1/2. The constant =fC(4i>) is obtained by evaluating both sides at Su — 
^fv + LUi. Thus one obtains 

A 2 = A = — (4 Su p(2v) + C(2 Su - 2v) + ((2 Su +2v)). (181) 
8a 

6.5 Centers of gravity 

Let z g2 denote the center of gravity of the area below the water-line and 4> 
denote the angle of z(u + Su) — z(u — Su) against the real axis in the complex 
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z-plane. Thus 

z{u+ 5u) - z(u~ du) = 2£e i,p . (182) 
An infinitesimal change of <f> yields 

2 

^2dz g 2 = - (z(u + 5u) — z(u — du))df, (183) 
which with d/ = £ 2 /2 and eq. (fT52"]) yields 



-4. 2 dz g 2 = 




(184) 



and thus 

z g 2 = =-r- (-ie 10 + const 2 ) . (185) 
0.A2 

Similarly one obtains for the center of gravity z g i of the area above the water-line 




(186) 



Thus both centers of gravity move on circles. This is a necessary and sufficient 
condition for the solution (see Auerbach[3] and references therein). Increasing 
u by 2a>3 will increase <f> by "0per ? and the center of gravity is rotated by f/'per- 
If it is not an integer multiple of 27r, then both constants vanish. If it is such 
a multiple, then continuity normally shows that it will vanish, too. We realize 
that then the line between the centers of gravity is perpendicular to the water- 
line and the distance h between both centers agrees with the expression (TTJ). 
This derivation did not use the solution of ([7]). 



7 Periodicity and Figures 
7.1 Limit of small e 

In the limit of small e one obtains (in this subsection we put 7*0 = 1) 

•m~^£±£ + <*">. - = ^ + o«»>. par, 

Thus for e = the limiting case of vanishing discriminant is approached with 
u>' tending to ioo and 

W3 = -^L= (188) 

from 

'_*_\ 2 2 



e ^^J-3- (189) 
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In order to determine oj' one may use the representation of the Weierstrass 
functions in appendix IA. 31 It is sufficient to use the terms up to linear order in 
q. With 

q = e 27ri ^i/^3 _ _ e mv'/u) 3 — —q 2 ^ q — e mu'/2u 3 (190) 

one obtains 

p{z) = (—^ + 8g (l-cos U )-±V u=— . (191) 

\ 2ui 3 J \ 1 — cos u 3 / 0J3 

For argument ^ + dz and small imaginary part of Sz one obtains 

,U>' „ , / 7T \ 2 / 1 , 7T(5z 



P ( T [-g-Sigsin— ). (192) 



One concludes from ei = p( " 2 " 3 ) 



9 = — 73; — " — 5T- (193) 
M(4 + M 2 ) 



Next v = ^- + 5v is determined from 



/'i = .^('')=-^T^-4 (191) 



12/i 
and yields 

sm( ) = - — . 195) 

u 3 2 

Now "0c, eq. (|115p can be evaluated with lo' = cu 3 + 2wi 

V> c = 47r + 8wiC(w 3 ) - 4iu> 3 C(2u) 

= 4?r + i(4w' + 8Sv)((lj 3 ) - 4iw 3 C(a/ + 2(5w) 

= 4?r + i(4w 3 + 8wi + 8<5u)C(w 3 ) - 4kj 3 C(w 3 + 2uj 1 + 2Sv) (196) 

= 4tt + 8i(wiC(^ 3 ) - w 3 C(wi)) + i(4w 3 + 8<5u)C(w 3 ) - 4iw 3 C(w 3 + 25v). 

The term 8i(wiC(w 3 ) — w 3 C(^i)) in the last line yields — 47r due to Legendre's 
relation (|334[) , the second one can be evaluated by means of eq. (|349[) , where q 
is now negligible and yields 

„ . / 7r ir5v\ . ( ttSv\ „ u 
ijj c = -27ricot - H = Tritan = ^ . - (197) 



x 2 w 3 / V ^3 / ^4 + /i 2 

Thus one obtains 

^ per = 2lu 3 = 2 7 r/i/v / 4 + M 2 = V'c (198) 

in the limit of small e. 
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Two dimensional bodies which can float in all directions are given by ip pcr = 
2n/n, thus for m = 1 and sufficiently small e. In this limit the Su can be 
determined from eq. (|14ip with 



<t(uj' + 5z) 



n , , „ . n ,nSz , „, , 

— {J +5z)-m-- tan( — ) + O(q), 



(199) 



2i 
— e 
nq 



-*» cos(^)e" 2 ^'+^) 2 / 24 + 0(g), (200) 



where eqs. p49|) and (|355[) and ^ — n , ^ — j2 nave been used. Then eq. 
(fT4Tj) yields 

tan(n 5u ) = n tan( Su ), (201) 

in agreement with the results obtained in refs. [7j,ll2}[8], where Su corresponds 
to £ — Sq and in ref. [6J, where Su corresponds to irp. 

A few cross-sections of the bodies are shown in figs. [10] to [23] For odd n the 
innermost envelope corresponds to density p = 1/2. 






Fig. 10 m/n = 1/3, 
e = 0.1 



Fig. 11 m/n = 1/3, 
e = 0.2 



Fig. 12 m/n = 1/3, 
e = 0.5 






Fig. 13 m/n = 1/4, 
e = 0.1 



Fig.* 14 

m/n = 1/4, e = 0.1 



Fig. 15 m/n = 1/4, 
e = 0.2 
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Fig. 16 m/n= 1/5, Fi S-* 17 Fi S-* 18 

e = o.l m/n = 1/5, e = 0.1 m/n = 1/5, e = 0.1 






Fig. 

m/n 




0.2 



Fig. 20 m/n = 1/6, Fi S-* 21 

e = 0.05 m / n = 1 / Q > e 



0.05 




Fig.* 22 

m/n = 1/6, e = 0.05 



Fig. 23 m/n = 1/7, 

6 = 0.1 



7.2 Periodicity 



In eq. (| 1 1 5|) an angle of periodicity ij} c has been denned. Here the periodicity is 
discussed for several regions in fig. 0J The angle of periodicity ^> per is defined 
as the change of the angle iji, as one moves from a point of extremal radius 
along the curve until a point of this extremal radius is reached again. Its sign 
is defined by the requirement that watching from the origin one starts moving 
counterclockwise. This yields 



A<0 



(202) 
(203) 
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Since due to eqs. P^lfT^l) 



dV> 
du 



c 



1 



one obtains 

-0 pcr = signr; Alp. 
Aip is obtained from eq. (j 1 09(1 by the increase of 

arger(u — 3u) — arg(j(w + v) + 2it3£(2v) 



(204) 



(205) 



(206) 



as u is increased by 2u>3. As shown in (|33H 1332)) this depends on the imaginary 
part of the argument of the functions a. Thus in addition to 



AV-o = 8viC(w 3 ) - 4iw 3 C(2w) 



(207) 



one obtains extra multiple of 2n. In region A' one obtains 3/T from o~(u — 3t>) and 
7r from er(w + u). This changes as one crosses to A", since there 3w < u//(3i). 
This can be seen from the sign change of p'(3v), eq. ([5^)1 due to the sign change 
of fj. One observes, that now the curve passes at its minimum radius on the 
other side of the origin. In the table (|208|) the ranges of v are listed in the 
various regions of fig. [Hand the corresponding ip peT are given: 



region inequality 


V'pcr ^pcr ^pcr ^pcr 


3i ^ < 2i 
A" < 3?j < 
B' < 3u < if 
B" if < 9fw < if 


V'c V'c 
V'c — 2vT 2/T — V'c 

V'c — 2vT 2/T — V'c 4/T — V'c 4/T — V'c 
V'c — 2vT 2/T — V'c 2/T — Vc V'c — 2/T 



(208) 



As examples figures out of the various regions are shown for Vper = 0. One 
of them looks like an eight. 



oo 




Fig. 24 m/n = 0/1, Fig. 25 m/n = 0/1, 

e = 0.5, e = 1.5, 

region A' region A" 



Fig. 26 m/n = 0/1, e = 2.5, 
region B' 
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Fig. 27 m/n = 0/1, 
e = 3.5, region B" 



Here and in the following the parametrization (|146p ip c — 2tt^ is used. It 
should be noted, that v, p(i>), and ip c vary continuously across all the borders 
between the regions A' , A", B' , and B" . ip c approaches for fixed e the limits 



Probably (I did not really check this) tp c is a monotonic decreasing function of 
e 2 for fixed e. 

7.3 More figures 

In the following more curves are shown. Now larger values of e are chosen 
as already done for the case m/n = 0/1. The next examples are for m/n — 
1/3, 1/4, 1/5, and 1/6. One observes that one obtains two branches for e > 5/2 
and m/n = 1/3 (figs. [33]-[33]). 




Fig. 28 m/n = 1/3, e = 0.9 Fig. 29 m/n = 1/3, e = 1.0 



lim ip c = 27r, 



lim ijj c = — oo, 

/j— >0,e<2 



lim ib c — — oo. 

e^e-2,£>2 



(209) 
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Fig. 30 m/n = 1/3, e = 1.1 Fig. 31 m/n = 1/3, e = 1.1 




Fig. 32 m/n = 1/3, e = 1.5 Fig. 33 m/n = 1/3, e = 2.6 




Fig. 34 m/n = 1/3, e = 3.0 Fig.* 35 m/n = 1/3, e = 3.0 
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36 
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Fig.* 53 m/n = 2/5, e = 0.1 Fig. 54 m/n = 2/5, e = 1.5 




Fig. 57 m/n = 4/5, e = 0.1 Fig. 58 m/n = 4/5, e = 1.5 



Note that the curve in fig. [55] obeys the convexity condition (|3"3")l (e 2 = 
—6.032657). There are no non-trivial solutions for the chord for m/n = 4/5. 
Some examples for negative m/n are: 
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Fig. 59 m/n = 

6 = 1.0 



-1/3, 




Fig. 60 m/n = -2/5, e = 1.0 




Fig. 62 m/n = -1/5, e = 0.2 




Fig. 64 m/n = -1/5, e = 0.5 




Fig. 61 m/n = -2/5, e = 1.0 




Fig. 63 m/n = -1/5, e = 0.2 




Fig. 65 m/n = -3/7, e = 0.4 
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Fig. 66 m/n = -1/7, e = 0.4 



Fig.* 67 m/n = -1/7, e = 0.4 




Fig.* 68 m/n = -1/7, e = 0.4 Fig.* 69 m/n = -1/7, e = 0.4 




Fig. 70 m/n = -1/2, e = 0.4 Fig. 71 m/n = -3/2, e = 0.4 

For m/n = 0/1 one obtains the Eight, which was already shown in figs. [24] 
to [571 Here some of the envelopes 7 are shown. 




Fig.* 72 m/n = 0/1, e = 0.5 Fig.* 73 m/n = 0/1, e = 0.5 
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Fig. 74 m/n = 0/1, e = 2.5 Fig. 75 m/n = 0/1, e = 2.5 




Fig. 78 m/n = 0/1, e = 2.5 Fig. 79 m/n = 0/1, e = 2.5 

Another example is shown in fig. [5U 

8 The Limit Case fi = 

The discriminant T> vanishes in the limit /i = and the elliptic functions degen- 
erate to trigonometric and exponential functions. In this limit two of the radii 
Ti coincide. 
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8.1 The case e > 2 

In this region one obtains with 

r4 = r (l + e), r 3 = r (l-e), r 2 = n. = -r (210) 
the expressions for e.; 

and the functions 

2A 2 

p(tt) = — + A 2 cot 2 (Au), (212) 
C(u) = — - + Acot(Au), (213) 

cr(u) = ie A2u2 / 6 sin(AM), (214) 
A 



A = *?E±. (215) 
2er V ; 



In this degenerate case one has 



From 



one concludes 



. = - = 7 _ ; w = oo. (216) 



( e + 2)( e +10) 
pM = 12AT~ (217) 



cot.(Ar) : cos(2Aw) = |, sin(2Aw) = i— ^— ^ (218) 



and 



= 2tt ( 2 - -^_= ) . (219) 



Evaluation of eq. (|109[) yields in a first step 



P r qju - 3v) c 2 C (2v)u 

2a 2 (2v) a(u + v) 



(220) 



r sm{\(u - 3v)) c2uXcot(2Xv) 
sin(A(w + v)) 

Multiplication of numerator and denominator by sin(A(u — v)) and use of sin a 
sin b = i(cos(a — b) — cos(a + bj) yields 

., _ _ cos(2AtQ - cos(A(2u - 4tQ) 2uX cot(2Xv) 

cos(2A W ) - cos(2Au) 6 ' ( ' 
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From eqs. (|218l) one obtains 



e 2 e^fe T ~^l 
cos(4Av) = — - 1, sin(4Av) = i— . (223) 



This yields 



= (e 2 - 2) cos(2Au) + ie^^I sin(2A^) - 6 /ro 

e-2cos(2Au) ' V ' 

Starting from eqs. (| 141 [) and (|154|) and using the expressions (|212[) to (|214[) 
one obtains the condition for the line segment of constant distance and its length 

tan(2Afe) = 2Ar tan( — - ^—^), (225) 

r 2 



'o 



l + 4^cot 2 (2Afe; 
Examples of the curves T are shown in figs. I80ltol83l 



(226) 




Fig.* 80 m/n = 1/3, e = 2.5 



Fig.* 81 m/n 

2.213594 



1/3, 
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Fig. 82 m/n = 1/4, e = 2.437087 



Fig.* 83 m/n = -2/1, e = 2.06559 



8.2 The case e < 2 

In the case of e < 2 one has 

r 4 =r (l + e), r 3 = r 2 = -r , n = r (l - e). (227) 
Then one obtains 

e 2 - 4 4 - e 2 



and the functions 



6e 2 r 2 ' u 12e 2 r 2 



o\2 

p(u) = — — + A 2 coth 2 (Au), (229) 

O 

C(«) = — — + Acoth(Au), (230) 

a{u) = i C - A2u2 / 6 sinh(AM), (231) 
A 



2 



A = (232) 
2er V 7 



In this degenerate case the half-periods approach 



a, 3 = oc, Wl = - = - ? _. (233) 



Again one obtains (|217[) and concludes 



coUkAc) = -i^/ 2 -^-!, cosh(2Au) = |, sinh(2A«) = 6 . (234) 



45 



Eq. p09p is evaluated in a first step 



2<7 2 (2v) a(u + v) y ' 

sinh(A(w + «)) y J 

Multiplication of numerator and denominator by sinh(A(it — v)) and use of 
sinhasinhb = i(cosh(a + b) — cosh(a — bj) yields 



cosh(A(2u - Av)) - cosh(2A«) r 2uX coth(2A „ 

i ; ' 

cosJ 

From eqs. (|234l) one obtains 



cosh(2Aw) - cosh(2Aw) ' 1 ' 



2 

This yields 



e 2 e^/? - 4 
cosh(4Av) = — — 1 , sinh(4Av) = i— . (238) 



= r (2 - e 2 ) cosh(2A^) + ieV4 - e 2 sinh(2Au) + se iu/ro 

° 2cosh(2Aw) - se ' 1 ' 

with s — + 1. Whereas for e > 2 the radius of the curve oscillates between | ?™3 1 
and f4, there are two curves for e < 2, one with largest distance = ro(l + e), 
the other with minimal distance |ri| = ro|e — 1| from the origin. This second 
curve is obtained by adding lj\ to v, or equivalently 7ri/2 to Xv. Thus cosh(4Aw), 
sinh(4Aw), coth(2Aw) are unchanged, but cosh(2Aw) changes sign. Starting from 
cq. (|237[) one finds that for the second curve one has to choose s = — 1 in eq. 
(|239p . Both these curves approach asymptotically the circle with radius r . 
The condition for the line segment of constant length and its length read 

tanh(2A(5u) = 2Ar tan(— - ^—^), (240) 

r 2 

If If 
0<- = . = < e, (241) 

r ° ^/(4-e 2 )coth 2 (2A(5u) + e 2 

and that for the line segment connecting both curves 

tanh(2A<5u) = — !— cot(— - ^—^), (242) 
2Ar r 2 

e<- = 26 < 2. (243) 

r ° ,J{4- e 2 )t&nh 2 (2\5u) + e 2 

Curves for e = 1.6 are shown in figs. [84] to [87] Figure [87] is shown in the 
movie 13] for a different solution 8u . 
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Fig. 84 m/n = 1/0, e = 1.6 



Fig.* 85 m/n = 1/0, e = 1.6 




8.3 The limit e = 2 

In the limit e = 2 one has 

r 4 = 3r , r 3 = r 2 = n = -r - (244) 

This yields 

ei = e 2 = e 3 = 0. (245) 
The Weierstrass-functions are degenerate to 

p(«) = -a, C(«) = - ff(u)=« (246) 
and there is no periodicity left, w = a/ = oo. With p(u) = — 1/Vq one obtains 

u = — (247) 
ro 

and 

^ = _ u-3iro e _ iu/ro = 3r + m e _ Wro _ 
w/r + i 1 - iu/ r o 

Again this curve approaches asymptotically the unit circle. The condition for 
the line segment of constant length reduces to 

i(-2Su/r +x-x) = / 24 gx 

l + i<5u/r ' V ; 
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cquivalcntly this may be written 

5u 



and I is given by 



tan 



e = 




5u 
ro 



X - X 



7q + 8u'' 



(250) 



(251) 




Fig.* 88 m/n = 1/0, e = 2 



Fig. 89 m/n = l/0 ; e = 2 



8.4 Elementary construction 

The circle 

I = r c 



-iu/r 



(252) 



is a second solution for e > 2. For e < 2 this circle is approached asymptotically 
for u — > ±oo. In all three cases the length of the line segment between this circle 
and the non-circular curves given above is 



21 = er . 



(253) 



This property allows an interesting construction of these curves. Consider 
a straight line (indicated in the figures by cyan and green) of length 2£ = ero ■ 
One end of it moves along the (magenta) circle of radius ro in figs. l90l to l96l 
The middle point of this line moves always in direction of the line itself, as if 
it were carried by a wheel aligned in this direction, creating the red envelope. 
Then the other end-point moves along the curve T. 

If e < 2, then there are two such curves, one outside and one inside the circle. 
Both approach asymptotically the circle for u — > ±oo. If e > 2, then there is 
only one curve, which for e = 2 still approaches asymptotically the circle with 
r — r oc 1/w 2 , whereas for e < 2 the difference r — r decays exponentially with 
u. For e > 2 there is only one curve, which oscillates between radii and r^, 
which allows for periodically closed curves. 

Examples for e > 2 are shown in figs. [90l to [93l below. 
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Fig.* 90 m/n = 1/3, e = 2.5 



Fig.* 91 m/n = 
2.213594 



-1/3, e = 




Fig. 92 m/n = 0/1, e = 2.309401 Fig.* 93 m/n = -2/1, e = 2.06559 

Examples for e < 2 are given in figs. l94l and l95l showing the outer and inner 
curve. 




Fig.* 94 m/n = 1/0, e = 1.6 Fig.* 95 m/n = 1/0, e = 1.6 

Finally the construction for e = 2 is shown. 




Fig.* 96 m/n = 1/0, e = 2 
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9 Linear Case 



In this section the limit is considered, in which tq tends to infinity, but the 
differences between the radii n is kept finite. In subtracting this average radius 
by taking the limit 

lim (z(u) - r ) (254) 

r ^oo 

one obtains curves, which repeat under translation in y-direction. Therefore 
this is called the linear case. One can distinguish three cases. 



9.1 Two branches 

The choice 

ta = r + d, r 3 = —tq - d, r 2 = -r + d, n = tq — d, d > d > (255) 
yields 

a = — , (256) 

2r (d 2 -d 2 ) 

d 2 +d 2 + 6dd 1 . . 

63 = 3(d 2 -d 2 ) 2 ~ 12rj(' (257) 



(P+d 2 - 6dd 1_ 

3{d 2 - d 2 ) 2 ~ 12r^ ' 



* 2 = —o ^ (258) 



-2(d 2 + d 2 ) 1 , . 

ei = — i + 259 

3(d 2 -d 2 ) 2 6rg' v > 



1 

v lies infinitesimal close to u>i/2, 



P( 2v ) = ^-^2- ( 26Q ) 



P'(f) = -^U + 0(1/^), (262) 
P(«) = ^f ) -^b) + ° (1/r ° ) - (263) 



This yields 



W =^+J V , 5t,=_i^— ^. (264) 
One starts from 

z = r 4 ^ffi e^W. (265) 
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Due to the factor in front which is of order rg , one has to take into account 
all terms of order 1 /tq in the fraction and in the exponential on the right hand 
side of this equation. With 

a(u-3v)=a(u-~uj 1 -3Sv) = -cr(u + ^uj 1 -3Sv)e 2 ^ UJl ~ u+35v '>^ 1 '> (266) 
and 

a[u + -wi - 3 Sv ) = a(u + v - 4 Sv ) = a{u + v)e~ i Sv ^ u+v ^ (267) 
one obtains 

z = r4 e 2tl (C(2«)-C(^i))+4<5i> (C(5Wi)-C(«+5Wi)) < (268) 

Further 

C(2v) - C(^i) = C(<*>i + 25v) - C(wi) = -2<5up(w 1 ) = -2 6vex (269) 

and 

(270) 



C(« -v)- C(« + w) + 2C(w) - p'(«) fa ( « 
yields 



z — r4e 4Sv (-eiu+C(^a)i)-C(Tt+5Wi)) (271) 
= r + d - 2r <5u f f (« - «) + CO + «) - ; t ^7 ; + 2e lU ) . (272) 

Adding a constant j/o to y, which corresponds to the angle \ with rox = yo the 
final result is 

2d(d 2 -d 2 ) , nrm . 

x-r = ±dT s — ^ , 273 

(d 2 -d 2 ) 2 p(u)-±(d 2 -5d 2 ) 

d2 ~d 2 f„ 1 . , 1 . 4(<?+«P) V 
» = yo + -2~ (C(« - + C(« + 2^0 - 3( , 2 _ j2) > j(?74) 

where the upper signs apply. 

Instead of starting from one can start from r±. Then the signs for the 
expression for x — r$ have to be changed, as indicated by the lower signs in 
equation (|273| . 

By inserting the extreme values of p(u) one finds that x — ro oscillates 
between d and d for one curve and between —d and — d for the other curve. 
This yields the two branches. Shifting the arc parameter by 2103 shifts the curve 
by y P er in the y-direction and leaves x invariant, 



x(u + 2uj 3 ) = x(u), y(u + 2uj 3 ) = y(u) + y por , (275) 

A(d 2 + d 2 ) 
3{d 2 -<P) 



A(fj 2 A-d 2 ) 

VP* = - J„ ? > 3 + 2(rf 2 ^d 2 )C(^ 3 ). (276) 
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In order to obtain the condition for the line segment of constant length one 
starts out again from eq. (| 141[) and expands both sides up to order O(l/ro) 



e i(x-*) = i + i ^L^£ ; (277) 

((2v)+((2v) = ((u>i + 26v) +((u) 1 + 2nw 1 +2 5v) 

= 2(n+l)C(wi)-4ei£v, (278) 

/ y-n ^i 2 ^ +V + V) _ e 2(l+n)(2 6u-2 5v)C{ui) 

cr(2 Su — v — v) 

a(2Su +(n+l)u 1 + 2 6v) 
a(2 5u + (n+ l)u! - 2Sv)' 



a(2Su + (n+l)u 1 + 2 6v) = ^c^+Cn+i)^) 
a{2 Su + (n + l)wi - 2 Sv ) 



(280) 



Then the terms in leading order cancel and the terms of order O(l/r ) yield the 
condition for Su 

yo-yo = ^Z, + _ d 2 Su ~ (<? - d ' 2 )(C(2 Su + (1 + n)wi) - (1 + n)C(wi)) (281) 

which for one and the same branch (even n) yields 
4C d 2 + d 2 ) 

Vo-yo= , 2 , 2 - - ^)(C(2 5w + wi) - C(wi)) (282) 
and for a line segment between the two different branches (odd n) 

y ° y° = oL a 5u - (d2 - d2)c(2 (283) 

3(gt — or) 

From eq. (|15T[) one learns that the length of the line segment of constant length 
for the same curve shifted by some yo — yo cannot exceed d — d. Examples are 
shown in figs. l97ltol99l 

These line segments are not the trivial ones parallel to the y-axis of length 
y por and integer multiples. The length of a line segment of constant length 
between two different curves has a lower bound d + d, but no upper bound. 
Indeed one obtains an infinity of solutions. A pair of curves is shown in figs. 
H00l to [T03l with three different lengths of the line segment. 

The curves shown here are also trajectories of massive charges in a perpen- 
dicular magnetic field. In this linear limit one obtains from eqs. E2J 12551 
1256)) that B increases linearly 

Amw , , 
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Fig. 97 d/d = Fig.* 98 d/d = Fig. 99 d/d = 
0.3 0.3 0.5 




Fig. 100 Fig.* 101 Fig.* 102 Fig.* 103 

d/d = 0.3 d/d = 0.3 d/d = 0.3 d/d = 0.3 



Note that on both trajectories in figs. 11001 to 11031 the charges move in the 
same direction. Such trajectories as well as those described in the following 
subsections have been obtained by Evers et al [2]. An example of a pair of 
drifting orbits is shown in fig. 1 of their paper. 

9.2 The limit d = 

In the limit d = one obtains 

e 2 = e 3 = (285) 
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e, = (286) 
p(u) = + ^coth 2 (^), (287) 

C(«) = -^ + ^oth(J, (288) 

*" ro = (289) 

2/ = y -" + dtanh(2u/d), (290) 
dni 

loi = — , w 3 = oo. (291) 

The curve is no longer periodic. The conditions for line segments of constant 
length are 

_ J — dtanh(2 Su/d), same sign of x — ro , , , 

m-Vo-^u-^ _ dcoth ( 2 5m y d j j different sign of a: - r . ( j 

Thus there is no chord of constant length for one and the same curve. The 
length of the line segments is given by 

2l=\y -y + 2Su\. (293) 

There are always two solutions for line segments between pairs of curves. 



Fig.*104d = Fig.* 105 d = Fig.* 106 d = Q Fig.* 107 d = 

A similar elementary construction as in subsect. 18.41 is possible. Now one 
end of the straight line of length 21 moves along the straight line x — ro = as 
shown in hg. 11071 Auerbach[3j has composed his examples for p — 1/2 from 
several copies of the central part of the curve and the corresponding piece of the 
straight line of fig. 11071 
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Evers et al. p3] show in the upper fig. 2 the variation of the classical 
trajectory of an electron as the gradient of B is varied, which corresponds here 
to a variation of d 2 /d 2 . The drifting orbit transforms into snake states, which 
will be considered in the next subsection. The limit case is the state d = 
shown here. 

9.3 One branch 

In the case of one branch (two fj are real, two are complex) the choice is 

r4 = r + d, r 3 = r - d, r 2 = -r + is, n = -r - is. (294) 

The calculation is rather similar to that in subsection 19. II Indices 1 and 3 have 
to be exchanged and d becomes is. This yields 



a = 5T) (295) 

2r (d 2 + s 2 ) y ' 

d 2 - s 2 + 6ids 1 
Cl ~ 3(d 2 + s 2 ) 2 ~12^ : ( ° 6) 

d 2 — s 2 - 6ids 1 
G2 = 3(d 2 + s 2 ) 2 ~Yh%' ( > 

2(s 2 ~d 2 ) 1 

« - + h (298) 

p{2v) = e 3 --V (299) 
Now v lies infinitesimal close to u'/2. One determines 

>j/ ^2 r j2 

P(y) = ^T^ + Oil/rl), (300) 

= - (<P ^, 2)a + ( 1 A3). (soi) 

P(«) = P^T 2 , +0(1Ao)- (302) 

2 ^o(« + s J 

From this one concludes 

w' d 2 + s 2 

V=^-+6v, 5v = -i^-ti-. (303) 

2 4r 

In the following calculation one has to replace u>\ by uj' and ei by e 3 and one 
obtains in analogy to (|272p 

z = r + d - 2r Q Sv (((u - v) + ((u + v) . f ^ . , + 2e 8 u) . (304) 

v PW ~ PW y 
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Fig. 108 s/d = Fig-* 109 Fig.* 110 

l§ s/d — 1.5 s/d =1.5 



Thus the final result is 

2d{d 2 + s 2 ) 

X ~ T ° ~ {d 2 + s 2 ) 2 p(u) + \{s 2 +hd 2 Y 

d 2 +s 2 A, 1 i\ „ 

y = vo + — ^ — [C(u - -u ) + ({u + 
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x — ro oscillates between d and —d. Moreover one has 
x(u + nw 3 )-r = {-) n {x(u) - r ), 



y(u + nuj 3 ) = y(u) + -y per , 



(307) 
(308) 



Vpc 



4(d 2 - s 2 ) 
3(d 2 + s 2 ) 



lo 3 + 2{d 2 + s 2 )CM- (309) 



The condition for the line segment of constant length can be obtained similarly 
as for the case of two curves. One obtains 



Va-Vo = - 3^ ~ 5 fa + (d 2 + S 2 )(C(2 fa 



(310) 



where one replaces C(2<5m +cj') — C( w ') = C(2fa +^ 3 ) — C( w 3), since and u> 3 
differ by a period. 

The examples of figs. 11081 to fTTTI show curves along which a bicyclist could 
drive. One and the same curve T for the front tire corresponds to different traces 
of the rear tire for bicycles of different lengths. 

In figs. U12l to [TTil the traces are more artistic, since periodically the bicyclist 
has to move back and forth with the rear tire, but always forward with the front 
tire. 






Fig. 112 s/d = 0.7 Fig. 113 s/d = 0.7 Fig. 114 s/d = 0.7 

If s/d decreases further the loops overlap stronger until at s/d — 0.458786 
the charge moves in an eight shown in fig. 11161 For this value of s/d one has 
J/pcr = 0. j/per becomes negative for smaller values of s/d. The expression (|284[) 
for the magnetic field reads now 

„ Amw , , 

B= -w^) {x - rQ) - (311) 
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These eights have still the property that between two of them one has line 
segments of constant length. When s/d becomes even smaller then the charges 
move into the opposite direction. In fig. 1 of ref. [H] two snake state trajectories 
are shown moving in opposite directions. In fig. 2 of this paper the eight can 
be seen where the electron is reflected in the magnetic bottle neck. 




Fig. 115 Fig.* 116 s/d = Fig. 117 Fig. 118 

s/d = 0.5 0.458786 s/d = 0.4219 s/d = 0.3 



10 The Carousel by Montejano, Oliveros and 
Bracho 

Oliveros, Montejano and Bracho[5l 0] have considered special cases of these 
curves. They consider a carousel, which is a dynamical equilateral (not equian- 
gular) polygon in which the midpoint of each edge travels parallel to it. Then 
the trace of the vertices constitute the curves T whereas the midpoints outline 
the envelopes 7. They perform explicitly calculations for pentagons. In general 
they allow for five different curves for the vertices. 

Denoting the curves by indices j one obtains solutions for the carousel from 
the curves introduced here by defining 

Zj (u) :=e 2: > iSx z(u), (312) 

and requiring 

\zj(u + Su ) — Zj-i(u — 5u )| = 11, z^(u + 10 Su) — zo(u). (313) 

The first condition is fulfilled with e 2i5x — r.h.s. of (|141[) . The second condition 
which guarantees that the pentagon is closed yields 

z 5 (u+ 10 (5m) = e 10iSx z (u+ 10 6u) = z a (u). (314) 

Thus 10 5u has to be an integer multiple of 2ct>3, 

k 

Su = —UJ3 (315) 
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with integer k. Then the condition reduces to 

e i(105 X +feV' pcr ) _ 1 



(316) 




An example is fig. 6 of ref. [5J. Since the curves go through the origin, 
one has e = 1. One realizes that the curves seem to be close to those for 
m/n = -1/3, e = 1 in fig. [5J which yields e 2 = -0.034862. One finds that 
condition (|316[) is fulfilled for e 2 = —0.025835 and apparently for two different 
k and two different lengths 2£ of the line segments. This is shown in figs. 11191 
and 11201 The second one seems to reproduce fig. 6 of ref. [3J. 

A solution for five 'Eights' is obtained for m/n — 0/1, e = 0.689454 in figs. 
11211 and 11221 The ratio of the arc length of the outer and the inner loop of the 
eights is 2/3. Fig. 8 in [5j shows the same eights (to the extend one can judge 
this with the eyes), but the pentagon shown there is misleading. It should be 
as in fig. EH 




Fig. 121 m/n = 0/1, Fig.* 122 m/n = 0/1, 

e = 0.689454 e = 0.689454 



The corresponding carousel for the linear case is shown in fig. 11231 and 11241 
This solution looks like that of fig. 2 of [4]. 
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Fig. 123 s 2 /d 2 = 0.0315653 



Fig.* 124 s 2 /d 2 = 0.0315653 



Of interest is the case where all five figures fall onto each other and the curve 
is closed. Then 

2tt k'nu 3 
iPper = — , ou = — - — (317) 
n 5 

with integer k' . One with n = 7-fold symmetry is obtained for e = 0.459104 in 
figs. 11251 11261 and should be compared with fig. 9 of [5 and fig. 3 of [4], the 
other one with n = 12-fold symmetry for e = 0.5727853 in figs. [T27l[l2iand is 
to be compared with fig. 1 of [4]. 
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Fig. 125 m/n = 1/7, Fig.* 126 ni/n = 1/7, 

e = 0.459104 e = 0.459104 




Acknowlegment I am indebted to Serge Tabachnikov for useful correspon- 
dence and to Peter Wolfle for useful discussions. 

A Weierstrass Functions 

A.l Definitions and relations of Weierstrass functions 

The Weierstrass function p is denned by 

p l2 (z) = 4p 3 (z) - g 2 p(z) - . 93 = 4(p(z) - e 1 ){p{z) - e 2 )(p(z) - e 3 ) (318) 
with the requirement that one of the singularities is at z = and 

ei+e 2 +e 3 = 0. (319) 
Commonly the two integrals are defined 

C(o) = - - / (P(z) ~ ^)dz, (320) 
a J z z 

a{a) = aexp (J^ (C(z) - ~)d^ , (321) 
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The function p is an even function of its argument, Q and a are odd functions. 
The Laurent and Taylor expansions start with 



p(«) = ^ + §« 2 + §« 4 + - 



C(«) = \- 9 i/-^~- (^3) 

*w - a -io a5 -io a7 -- ^ 

There exist addition theorems 

«. + ») - - pW -p W + gW^(|! , (325) 

«. + ») = CW + CW + ^''°i" p, ff , (326) 

2 p(a) - p(&) 

<r(a + b)a(a - b) = -a 2 {a)a 2 {b){p{a) - p(b)). (327) 

If to is a half-period, that is u> itself is not a period of p, but 2lu is, then the 
following relations hold for integer n 

p(a + 2nu>) = p(a), (328) 
({a + 2nu) = C(a) + 2nC(u>), (329) 
a{a + 2ncj) = {-) n <j{a)e 2n{a+nuj)a ^ . (330) 

The increase of the argument of c(a), if the real part of a is increased by 2u>3, 
is obtained from eq. (|330[) 

argcr(a + 2^3) — arg a(a) = — 7r(2m + 1) + 23aC(^3) (331) 

with integer to, where — ir(2m + 1) comes from the minus-sign in (|330[) and to 
has to be determined. If a is real, then a changes sign, whenever a is an integer 
multiple of 2^3. If a has a small positive (negative) imaginary part, then a passes 
the zero clockwise (counter-clockwise) in the complex plane, yielding a decrease 
(an increase) of the argument by ir. Further jumps by 2tt occur, whenever the 
imaginary part of a crosses an integer multiple of 3wi . Thus to in eq. (|33ip is 
given by 

2to3wi < 3a < 2(to + (332) 

From (|326[) one obtains 

C(« + b) + C( a -6) = 2C(a)+ ( P ' {a) (hV (333) 

p(a) - p{b) 

Legendre's relation reads 

rjuj' - Wu = — for 3— > 0, (334) 
2 u> 

where 2lo and 2co' span one elementary cell. 
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A. 2 Double-periodic functions 

Two basic theorems on double-periodic functions are 



Theorem I Two elliptic functions which have the same period, the same poles, 
and the same principal parts at each pole differ by a constant. 



Theorem II The quotient of two elliptic functions whose periods, poles, and 
zeroes (and multiplicities of poles and zeroes) are the same, is a constant. 
Due to this second theorem the elliptic function 



f 



U 1 ,U 2 —U m 
Vi,V 2 ...V n 



(335) 



is defined as the periodic function with zeroes at Hi, u 2 , ■■■ u m and poles at 
vi, V2, ■■■ v n in an elementary cell. Zeroes and poles at the origin are counted 
separately. If n > m, then the function has an n — m-fold zero at u = 0. If 
n < m, then the function has instead a pole of order m — n at u = 0. The 
function £ is normalized so that 



£ n - m {u) = u n ~ m + 0(u 



n— m+1 



(336) 



Moreover one has to require that the sums Ui and Vi differ only by a 
period 2w 



i=l 



i=l 



2 



(337) 



with integer kj. The functions £ can explicitly be expressed in terms of the 
functions a 

Ui,U2—U r , 
Vl,V2--V„ 



(338) 



Using eqs. p30l) and p34p one finds that this function is double-periodic. g 2 
and #3 are not explicitly indicated. In these eqs. it is always assumed that they 
are the same in all functions. 
Multiplication yields 



£n+n' —m~m' I ^ 



Ul,U 2 ...U„ 

Vi,V2—V n 

Ui,U 2 



Im-, u 2' 



v ll v 2 ...v n ,v' 1 ,v 2 . 



u 1 ,u 2 
v[,v' 2 

:U'„ 



■v' 



If locations of zeroes and poles coincide they have to be cancelled 



f 



oi, a 2 , 
ai,a 2 , 



., a k ,ui,u 2 -.-u m 
.., a fc , vi,V2—v n 



Ui,U 2 ...U m 

v 1 ,v 2 ...v n 



(339) 
(340) 

(341) 
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Examples are 



p(u) - p(v) 

p(u) - p(v 3 ) 
p(u) - p{v) 



£-2 yu 

£ I u 



V3, -V3 

v, —v 



a(2u) 
in)' 



(342) 
(343) 
(344) 



and 



1 p{u) p'(u) 
1 p(v) p'(v) 
1 p(w) p'(w) 



= -2E-n u 



v, w,—v — w 



£-2 V 



w, —w 



—2o(u — v)a(u — w)a(v — w)a(u + v + w) 
<7 3 (u)<t 3 (v)<j 3 (w) 

from which one concludes that the determinant vanishes, ifu + v + w = 



(345) 



1 p(u) p'{u) 
1 p{v) p'(v) 
1 p(u + v) -p'(u + v) 



0. 



(346) 



In the limit u — > v one obtains 



1 p(v) p'(v) 

p'(v) p"(v) 

1 p(2v) -p'(2v) 



0. 



(347) 



A. 3 Representations of Weierstrass functions 

The Weierstrass function ((z) with elementary periods 2lo and 2a/ has simple 
poles with residues 1 at 

z = 2muj + 2nJ (348) 

for integer m and n. It reads 

/ 00 / / \ 

C (z) = ^ + _( co t(- ) + ^(cot(- + - r)+ c t ( --- r) ) 



n=l 



-z+— sinu h4> — 

u 2lu \ 1 — cos u ^— ' 1 — 2(7™ cos 

\ n=l 



cos u + q 



2)i 



with 



7TZ 



U = 



q = exp(2i7ro/ / w) , 



(349) 



(350) 
(351) 
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\ n— 1 v cj / , 

= ^(f) 2 fl-8V , g " g V (352) 
The derivative of the ^-function yields the Weierstrass p-function 

+ (ir?(t) + § ( sin 2 (f + 2^1) + sin 2 (f - ^) } ) 

^ 2 /tt\ 2 / 1 g"(cosu-2 g " + g 2n cosM) \ 

w \2uJ il-cos?i ^ (1 - 2q n cosu + q 2n ) 2 h ' 

The constant term to p(z) fixes the constant term in the Laurent expansion of 
p(z) = 77 + 0z° + ... to vanish. The derivative of the Weierstrass function reads 

3 . / 1 



AZ) = - 4 (£) 3sinU ( 



(1 — cosm) 2 

4 9™ (1 + 2<?" cos u - 6<7 2 ™ + 2g 3 ™ cos u + q in ) \ 
2-^ (1 - 2(7™ cosu + 9 2 ™) 3 J ' 



71=1 

Finally the integral of £ is obtained 

<r(z) = zexp^(C(z')-3M" 



2w . u ( r\z 2 \ 1 - 2q n cos u + q 2n 
= V Sm 2 CX H^J n l-2 g » + g2 " • (355) 

v 7 n— 1 

These representations of the Weierstrass functions have been used for the cal- 
culation of the figures. 

B Symmetric Polynomials 

Four symmetric polynomials P q , P r , P m and P are introduced. Symmetric 
means that they are invariant under permutation of the Ti and qi, resp. Three 
of them are defined by 

P q = ^(<?i + 92 - 93 - <74)(<7i - 92 + 93 - 94)(9i - 92 - 93 + 94), (356) 

P r = ^(n + r 2 - r 3 - r 4 )(ri - r 2 + r 3 - r 4 )(n - r 2 - r 3 + r 4 ) 

= -(ri + r 2 )(ri + r 3 )(r 2 + r 3 ) = (r 4 + ri)(r 4 + r 2 )(r 4 + r 3 ), (357) 

Pm = ^*9j ■ = ^(91 + 92 + 93 + <a? -\{ql + ql + ql + ql)- (358) 
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The identity for P r is derived from J+ — 0. Expanding the last expression in 
(|357|) as polynomial and subtracting r 4 (r 4 + T\ + r 2 +7*3), which vanishes one 
obtains 

Pr = rir 2 r 3 + r x r 2 r± + rir 3 r 4 + r 2 r 3 r 4 . (359) 

Realizing that 

91 + 92 - 93 - q± = r\ + r 2 2 - r 2 - (n + r 2 + r 3 j 2 = -2(r 3 + ri)(r 3 + r a ) (360) 
and similarly for permutations one finds 

P 2 r = Pg. (361) 

Consider 

r 4 P r = r 4 (r 4 + ri)(r 4 + r 2 )(r 4 + r 3 ) 

= n(rf + r|(n + r 2 + r 3 ) + r 4 (r!r 2 + fir 3 + r 2 r 3 ) + rir 2 r 3 ) 

= 7^94(94 - 9i - 92 - 93) + rir 2 r 3 r4, (362) 
where use is made of 

r 1 r 2 + r 1 r 3 + r 2 r 3 = i ((n +r 2 + r 3 ) 2 -rf -r 2 -r§) = -(94-91-92-93)- (363) 
Using also permutations one obtains 

= (n + r 2 + r 3 + r 4 )P r = A ri r 2 r 3 r A - 4P (364) 

with 

^ = -^(9? + 9 2 + 9 3 2 + 9 4 ) + ^m 

- ^(9i + 92 + 93 + 94) 2 - J(9? + 9 2 +9,3+94)- (365) 

This yields 

r\r 2 r 3 r± = P (366) 

in terms of the squares qi = rf. 
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